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Abstract 

' We introduce and implement a new, extended approach to placing bounds on trilinear R- 

parity violating ) couplings. We focus on a limited set of leptonic and semi-leptonic processes 
^ involving neutrinos, combining multidimensional plotting and cross-checking constraints from 

CIh' different experiments. This allows us to explore new regions of parameter space and to relax a 

P> i' number of bounds given in the literature. We look for qualitatively different results compared to 

r~| I those obtained previously using the assumption that a single coupling dominates the R-parity 

violating contributions to a process (SCD). By combining results from several experiments, we 
' identify regions in parameter space where two or more parameters approach their maximally 

OO I allowed values. In the same vein, we show a circumstance where consistency between indepen- 

dent bounds on the same combinations of trilinear coupling parameters implies mass constraints 
among slepton or squark masses. Though our new bounds are in most cases weaker than the 
I SCD bounds, the largest deviations wc find on individual parameters are factors of two, thus 

^\ I indicating that a conservative, order of magnitude bound on an individual coupling is reliably 



o 



estimated by making the SCD assumption. 



X 

H : 1 Introduction 



It is well known that the standard model (SM) admits some "accidental" symmetries such as the 
separate conservation of baryon (B) and lepton (L) number. In other words, the requirement of 
gauge-invariance and renormalizability of the operators that appear in the Lagrangian does not 
allow the presence of terms that violate baryon or lepton number conservation. In the framework 
of the Minimal Supersymmetric Standard Model (MSSM) this is no longer true. In this model, 
operators that carry the same baryon and lepton number of the standard model, but different spin 
or mass dimension (the superpartners), violate B and L conservation, which can be enforced by 
hand with the introduction of R-parity. This additional discrete symmetry of the spinorial charges 
allows the lightest supersymmetric particle (LSP) to remain stable and is defined as [1]: 
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i?=(-l)3{S-^)-25^ (1) 

with S being the spin quantum number. All the standard model particles have R = 1, while their 
superpartners have R = —1. 

The phenomenological signatures of an unstable LSP have been investigated extensively in a 
variety of papers, both at lepton [2, 3] and hadron [4] colliders. Generally the signatures are the 
consequences of new interaction terms that arise in the superpotential or in the soft supersymmctry 
(SUSY) breaking part of the Lagrangian when the assumption of R-parity is lifted. Wide attention 
has been given to extracting bounds on these new couplings from precision tests of the standard 
model and from cosmological constraints. The extent of the literature on the subject is daunting: 
we refer the reader to Ref. [5] and references therein for a comprehensive review. 

Given the impact that a precise determination of the coupling size has on the phenomenological 
consequences, wc think it is important to obtain them in the greatest possible generality. In this 
paper, we do so by relaxing some of the assumptions that are commonly used in the literature. After 
reviewing the form of the R-breaking couplings and deriving the effective Lagrangians of interest, 
in Sec. 2 wc describe the assumptions commonly used in the literature and introduce our extended 
approach. In Sec. 3 we derive new bounds on the R-breaking couplings from leptonic processes, 
while the bounds from semi-leptonic processes are treated in Sec 4. In Sec. 5 we summarize our 
results and conclusions. We use new PDG2008 [6] data to obtain bounds at 2a that are, in cases 
where new data has become available, more stringent than the existing ones under the standard 
SCD assumptions. 



2 R-parity violating couplings and low energy effective Lagrangians 

There is no theoretical argument that prevents the superpotential from having the following bilinear 
or trilinear terms^: 

/t = Kjk^abLfL^jEk + Kjk^abLiQ'jD^ + X'lji^eimnUf''Dj^"^D^"' (2) 

and 

fs = f^'i^abLfH^, (3) 

where the carets label the superfields corresponding to the standard model fields, the indices i, j, k = 
1,2,3 label the fermionic generations, a,b = 1,2 are S'?7(2)-doublet indices, while l,m,n = 1,2,3 
are >S'[/(3)-triplet indices. The Ay^ couplings are antisymmetric in i,j due to the antisymmetry in 
a, b, imposed by SU{2), while the A^^j. are antisymmetric in j, k due to the complete antisymmetry 

^Here and throughout we follow the conventions and notations of Ref. [7]. 
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of Qmn; required by SU{3). One can see that the first and second terms in Eq. (2), and Eq. (3) 
violate L conservation, while the third term in Eq. (2) violates B conservation. On the other 
hand, phenomenological considerations show that the trilinear terms in Ajj^ and A^^^ cannot be 
simultaneously present with values large enough to affect the processes we study here, otherwise 
squark-exchange would lead to unacceptable rates for proton decay [8, 9]. 

Eqs. (2) and (3) introduce 48 new complex parameters to the MSSM: 3 dimensionful parameters 
from the bilinear couplings, and 9 + 27 + 9 = 45 dimensionless parameters from the trilinear. 
Along with the superpotential terms, B and L can also be violated by 51 additional soft SUSY- 
breaking terms in the Lagrangian. Since they are not pertinent to the following discussion we 
will not write them explicitly here. They can be found in Ref. [5], along with a discussion of the 
choice of bases in which the bilinear term in the ^ superpotential, Eq. (3) , is rotated away by an 
SU (4) transformation, so that the sneutrinos acquire a vacuum expectation value under electroweak 
symmetry breaking [10, 11, 12]. 

Consistent with the existing literature on trilinear ^ bounds as reviewed in [5] , here we choose 
to work in the mass basis, assume all bilinear ^ terms in the tree-level Lagrangian are absent and 
base our analysis solely on the trilinear terms. 

Since R-parity violating terms are neither forbidden by gauge invariance nor by renormalizabil- 
ity, but rather depend on phenomenological consistency, one can wonder to what extent R-parity 
could be broken, i.e. how big are the couplings appearing in Eqs. (2) and (3). Restricting our 
discussion to Eq. (2), determinations of the couplings' size are generally obtained in the literature 
by comparing an effective Lagrangian expressed in terms of the A, A' and A" couplings with the 
neutral and charged current interaction effective Lagrangian that describes fundamental tests of the 
standard model. We largely confine ourselves in this paper to flavor-conserving cases, to keep the 
presentation focused. The most general effective Lagrangian for fermion-fermion neutral current 
interactions II —>■ ff at low energies reads: 

jC = -4V2Gf + e'L)L + (5^ + e^)i?]//7^[(5{ + e{)L + {g{, + e{)R]f, (4) 

where Gp = V^g^ / {SM"^) is the Fermi coupling constant, L = (1 — 75)72 and R = (1 + 75)72 
are the chiral projectors, ql and are the coupling to the chiral components of the fundamental 
spinors, and the e's describe the "non-standard" part of the interactions. One requires that the 
R-breaking contributions do not exceed the limit imposed by the precision of the experimental 
measurements, thus obtaining bounds on the couplings. 

As we have mentioned above, the simultaneous presence of leptonic and hadronic R-parity 
violating couplings is tightly constrained experimentally by the stability of the proton. One may 
then choose to consider either Ajj^ couplings, or A^'-^ couplings to be negligible. In this paper we 
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deal strictly with processes that involve A and A', as their corresponding experimental signatures 
are clearer and the reported uncertainties are smaller. 

An effective four-fermion Lagrangian, applicable to processes at energies small compared to the 
weak scale, can be obtained from the superpotential of Eq. (2): 
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r.s 



dSrdSc: 



s=s 



5=5 



(5) 



where r, s span the superfields S of the superpotential, and 'tpr,s are the Majorana fermion fields 
entering the supermultiplets. The part involving semi-leptonic interactions is given by the second 
term in Eq. (2): 

(6) 



/ 3 >!i.j,{vidjD^ - liU.iD^ 



Application of Eq. (5) to this term yields 



+dLj1pUiL'lljjjC - ULji^eiL'tpi^c) + h.c. 



(7) 



The vertices can be obtained by defining Dirac spinors as 



(8) 
(9) 



so that one gets for the interaction part of the Lagrangian, 



->'ijk[dRk''f^dj + dLjdkLvi + DidkLdj - d^j^,^ef'Luj 
-eudkLuj - ULjdkLci] - X'*jf.[dRkdjRu[' + d\.ViRdk 
+u\djRdk - dRkUjRef - e'j^-UjRdk - n^^eii^djt]. 



(10) 



The effective Lagrangian for scalar mediated four-fermion interactions can be obtained by com- 
bining the vertices of Eq. (10) and applying Fierz identities to the result: 



1 _ - 1 _ 
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(11) 



The effective Lagrangian of Eq. (11) introduces 135 independent parameters: 9 combinations in 
any two of the indices times 3 combinations in the remaining index which runs through the families 
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of 5 possible exchanged sparticles. The leptonic interaction effective Lagrangian is obtained by 
applying the same procedure to the first term in Eq. (2): 



/ 3 Xijk{i>iejE^' - eiUjE'^). (12) 



One gets [3] 



C 
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eff - 2 

1 



1 _ _ 1 _ 

-^^{uat^Lvi){ej^^Lej) i^{eij^'Lui){i'jji^Lej) 

mi mi 

, ^Rk <^Rk 

{ej'j'^Lej){ekj^Rek) \—(pj-i^'Lvj)(ek-l^Rek) ) + ^ j) 



(13) 



where i < j is understood in Eq. (13). The same antisymmetry in the i and j indices of the A 
couplings reduces the number of effective independent couplings encompassed in Eq. (13) with 
respect to the semileptonic case. There are 45: 3 free parameters in a 3x3 antisymmetric matrix 
multiply 3 possibilities for the remaining free index that carries the dependence on the mij^^; 
plus, there are 6x3 possibilities with left-handed sneutrino exchange and 18 more possibilities with 
left-handed selectron exchange. 

The limits in the literature are obtained under the assumption that a single coupling dominates 
the R-parity violating contributions to a process (SCD). This assumption rests on the premise 
that some hierarchy exists between the leptonic, semileptonic and hadronic couplings, or between 
different fermionic families. Besides, the couplings often enter as sums of squares, so that one might 
guess the most conservative bounds follow from this hypothesis. We found that in most cases this 
is not so. 

It is an open question whether such a hierarchy does indeed exist. In the absence of a theo- 
retical guide, we apply a "multi-parameter" approach to placing bounds, to explore new regions of 
parameter space. In Sections 3 and 4 we give examples of our approach and contrast the results to 
those of the SCD simplification. 

2.1 Notation and conventions 

We think it is important at this point to clarify our notation, as the originality of our contribution 
rests in making explicit use of some properties of ^ -couplings that are often overlooked in the 
literature, partly because the established notation bears some elements of ambiguity. As far as 
SCD is concerned, the concept was originally formulated by Dimopoulos and Hall [2]. As generally 
applied, one assumes that a single coupling (or a single product of couplings) is much larger than 
the others which, therefore, can be neglected when placing bounds. As is the case in most of the 
corrections to the SM that involve R-parity violating couplings, more than one coupling is present 
and often this simultaneous presence is not clear in the notation. For example, when the process 
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at hand involves the four-fermion interactions described by Eqs. (11) and (13), the initial and final 
states of the scattering or decay are supposed to be completely known, whereas the exchanged 
sparticle, whether a squark or a slepton, can be of any generation. Thus, since this flavor is 
unknown, one has always to sum over the families of the sparticles compatible with the relevant 
vertices. So, it is important to understand that a bound that reads, for example, |Ai2fc| < 0.15{eRk) 
can be taken to mean either 

|Ai2fe| < 0.15 (y^^) , (strong version) (14) 

for each A; = 1, 2, 3, or 



(~^) - O-^Ol^ GeV-^ (weak version) (15) 

As in Eqs. (14) and (15), we adopt the standard 100 GeV scaling of sfermion masses throughout. 
One of the ways of implementing the SCD convention consists in setting all but one \i2k in Eq. (15) 
to zero, thus effectively obtaining Eq. (14). The strong version produces bounds that are obviously 
more conservative than the weak one, so we will display the form (14) every time we place a new 
bound on a coupling, with the caveat that the reader can interpret it in the form of Eq. (15). We 
will state explicitly when we make an exception to this rule. 

If the initial and final states of the four-fermion process involve the same vertices, the R- 
brcaking couplings enter the process only through their modulus squared. In the literature it is 
then customary to express the corrections to the SM as functions of simplified quantities: rjjfc(/j) 
(but also rijkilj), r,jk(lk}) or r[-^{fi) (but also r[-j^{fj), r[-f^{fk))- In light of what we have explained 
above, we want to make clear that these are symbols that stand in full for: 

where the scaling factor 4\/2Gf comes from the general form, Eq. (4). Thus, they admit a sum 
over the flavors of the exchanged sfermion /, (or fj, fk) which, depending on the case, can be a 
slepton (li) or a squark (gj). It is also clear that the value of the mass of the exchanged sparticle is 
always left unknown. If, instead, the SUSY process involves different vertices, then the correction 
to the standard model is expressed as a function of a product of couplings, of the kind Xijk ■ Ksk 
(equivalently, A^^^ • A^^^^ or X^jk ■ A^^^). In these cases too, a sum over k needs to be considered. 
Analyses that return one product as dominant are a common extension of the SCD. 

We have decided to label the fermion (sfermion) generations by a number index i (or j or k) 
= 1,2,3 whenever the families are summed over, as in Eqs. (15) and (16), or when one of the 
indices is free to take any values, as in Eq. (14). But, for clarity's sake, if the bound involves just 
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one single particular coupling we will label the generation by name so that, for example, ejn cr, 
^L2 t'/iL, ULz ^ iL, and so on. 

We have mentioned above that implementation (14) of the SCD produces bounds that are more 
conservative. In most cases, though, a physical process cannot be expressed in terms of only one 
combination of couplings such as (16). The bounds from experiment are placed generally on a 
function of combinations 

F (rijkik), rrstiL), r'l^niQn), -) • (17) 

It is a common approach in the literature to set all the r's of Eq. (17) but one to zero, so as 
to place bounds on the surviving combination of couplings. Moreover, one term at a time of the 
combination is then assumed to dominate. It is this particular implementation of the SCD that 
wc find excessively severe, as it reduces the dimensionality of the allowed regions of parameter 
space thus missing any information on the combined action of different couplings embedded in the 
function F. In the next two sections we show that allowing the full dependence on F does indeed 
give more information and in some cases also extends the allowed bounds on the couplings. 



3 Leptonic case 

In order to show how our approach works, we start with some classical examples in the leptonic 
case [3, 5]. We begin with constraints required by universality in muon and tau decays, then take 
up the constraints from v^e, fgC, and i^ee elastic scattering cross section measurements. 

3.1 Muon and tau decays 

Let us consider the two following ratios: 

and 

r(r IJ. U^Vr) 

which are sensitive to violation of lepton universality. By comparing the tree-level effective La- 
grangian of Eq. (13) with the SM, one can derive bounds on some of the A-couplings. The SUSY 
processes that contribute to the decay (18) are shown in Figure 1. Those for (19) can be obtained 
by replacing j = 2 ^ 3 in Figure lb. 

Besides, one also needs to consider the A-dependencc of the Fermi coupling constant Gp [2]. As 
is well known, Gj? is experimentally determined from measurements of the muon lifetime. Therefore, 
when dealing with R-parity violating SUSY, Gp receives a correction from the SUSY processes that 
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contribute to ^-decay, Figure lb. The correction is given by [3]: 

where a sum over repeated indices is intended, as explained in Sec. 2.1. We have also used the 
notation introduced in Eq. (16). 

Taking into account the processes of Fig. 1, Eqs. (18) and (19) can be expressed in terms of 
their SM expressions and yield, to first order in R-breaking couplings [3], 

R^^ = Rfl^' {1 + 2 [r23k{eRk) - rMeRk)]} (21) 

and 

Rr = Rf^ {1 + 2 [ri3fc(eRfe) - r23k{eRk)]} , (22) 

again, with the conventions of Eq. (16). As explained in the discussion preceding and following 
Eq. (17), if we were to use the SCD at this point, we would consider one r-combination at a time 
and obtain a bound on each of them when the remaining couplings are put to zero. By using 
the measured values of Rrfj, and Rr [6] and the standard model values after radiative corrections 
([13, 14] and References therein) for -Rf^^ and we would obtain at 2a: \X23k\ < 0.063 (eRk) 

and |Ai2fc| < 0.045 {cRk) from i?^^, |A23fc| < 0.051 (cRk) and |Ai3A;| < 0.048 (eRk) from Rr, where 
the dominant uncertainty is the one on the r lifetime, and we have used the conventions introduced 
in the discussion preceding Eq. (14). In principle, though reasonably well motivated and not 
inconsistent, there is no theoretical justification for considering one coupling at a time, or a sum 
over families at a time. As Figure 2 shows, the full dependence on the couplings presents a richer 
structure. 
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Figure 2: a) |A23fc| and |Ai2fc| are underconstrained at 2a if one takes solely into account 
Rr^,, Eq. (21). b) 2a bound region on X23k{eRk), \i2k{eRk) and Xi^kieRk) from i?^^, Rr 
and muon lifetime combined data. The allowed region can be enclosed in a box of size 
{|A23fc|(ei?fc),|Am|(ejjfc),|Ai3fc|(ejjfc)} < {0.075,0.043,0.082}. 



Figure 2a shows that Eq. (21) admits degeneracies on the couplings. When taken together, 
|A23fc| and |Ai2A:| can be taken arbitrarily large^, since they cancel each other. A similar picture 
holds for Rt, as Eq. (22) has the same form as Eq. (21). Thus, an approach extended beyond SCD 
consists in trying to limit and reduce those degeneracies by combining different experiments that 
involve the same couplings. 

In particular, the measurement of the muon lifetime can be used to determine a first bound 
on the sum of couplings Ai2fc, when a right-handed charged slepton is exchanged, Eq. (20). The 
result is dependent on radiative corrections and on the renormalization scheme. Expressions for 



Ai2fc/"T'g^j. can be derived in the on shelly MS and Novikov-Okun-Vysotsky (NOV) renormalization 
schemes [15], and are given by 



Al2A:P SGp 



M| \/2G'fp| scheme (sin^ Ow COS^ 6*^/ ) I scheme (1 - Ar [scheme) _ ^ 



7ra 



(23) 



where the scheme dependence of p, s\v? Ow and Ar can be found in Table 1. According to the 
different renormalization scheme, we find the lu-bounds on A12A: obtained from the muon lifetime: 

(on shell) \\i2k\ is excluded (24) 
(NOV) |A,..| < 0.011 (^) . (20) 



Consistent with the underlying perturbation expansion. 
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Scheme 


sin^ 9w 


Ar 


P 


On shell (o.s.) 




1 a Pt 

a(Mz) tan^ewlo.s. 


1 


MS 




1 - + -[16] 


1.01023(22) 


NOV 


1 (l 7Ta{Mz) 

2 \,4 ^GfM^J 


1 

a(Mz) 


1 



Table 1: Analytic expressions for sin^ 9w, Ar and p in the on shell, MS and NOV renormalization 
schemes, pt = SGprnf / 8 V2Tr'^. We use the following average values at la [6]: sin^ 9w\-ms ~ 
0.23119(14); a-^{Mz) = 127.925(16); a'^Mz) = 128.91(2); sin^^iylNOV = 0.23108(5); Ar|o.s. = 
0.0369(14); 

^^\~MS — 0.06962(12). The ellipsis indicates non-leading order terms that can be found 

in [16]. 
And at 2a: 

lOOGiv) ^^'^^ 

(NOV) \Xnk\ < 0.015 {j^^) . (29) 

We can therefore decide to use this bound to limit the degeneracies present in i?^-^ and Rr- Figure 
2b shows the 2(7-allowed region of parameter space in A23fc, Xi2k and Xi^k when the PDG2008 data 
for Rth, Rt and muon lifetime are combined. We use the MS bound, Eq. (28), on \i2k- One can see 
that the A-parameters undergo a extension up to a factor of two with respect to the value obtained 
using the SCD. When scaled to the masses of the exchanged sleptons the 2a region shown in Fig. 
2b can be enclosed in a box of size {\\23k\{eRk),\>^i2k\{eRk),\>^i3k\{eRk)} < {0.075,0.043,0.082}. 
The 2a bounds are: 

and 

|Ai3fc| <0.071 ""^^^ ]. (31) 

The other schemes yield similar values. Note that the new bound on Ai3jt docs not include A133, 
since this coupling is separately and more severely bounded by the mass [2]. Similarly, there 
exist strong bounds on many pair-wise products of the couplings above, coming from experimental 
bounds on decays disallowed in the SM [5, 17]. However, there are always combinations of Ai2fc, 
A13J; and A23jk that are still unconstrained by the bounds on products. This comment applies to all 
the cases that wc are considering. To list the detailed conditions takes us beyond the aim of this 
paper, so we leave them as implicit. 
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3.2 Neutrino - electron scattering 

Wc now turn to the flavor diagonal neutrino - electron scattering processes + e ^ f ^ + e and 
z^e + e ^ t'e + e. In the ly^e and z/ge examples, the energy is always large enough to neglect the 
electron mass in the kinematics, while in the case, the neutrino energies are in the MeV range, 
which requires us to keep the electron mass cflrects in the kinematics. The individual left -and right- 
handed couplings or, equivalently, axial and vector couplings, have been extracted individually in 
the experiments on the f^e case [18], making the analysis of bounds on R-parity violating parameters 
quite straightforward; we begin with this process. 

3.2.1 i/^ + e ^ z/^ + e 

Neglecting the terms proportional to the electron mass, the total cross sections for f ^ + e — f ^ + e 
and 17^ + e — ^ + e can be written as 

= ^ (^i + lal) (32) 

and 

a(l7,e) = ^(ffU^ffi)- (33) 

The direct-channel Mandelstam variable s = 2meEy in the target-electron rest frame. We can write 
gi and gn in terms of the weak angle and the R-parity violating parameters (Figure 3) as 

gL = gl'' -{l+gl'^yMeRk) (34) 

and 

9B. = Qr'^ + ri2i(eL) + r23i(fL) - g^^ ri2k{eRj,), (35) 

where xw = sin^ Ow, and gf^^ = xw — \ and g^ = xw are the SM expressions for the L and R 
couplings. 

Since the experimental averages, with errors, are reported by the Particle Data Group [6] for 
gA = gh — gR and gv = gL-\- gR, we use these forms to obtain the bounds on the R-parity violating 
couplings: 

SA = gl^ [1 - ri2kieRk)] - ri2i{eL) - r23i(fL) - r^kieRk) = -0.507 ± 0.014 (36) 

and 

gv = g^^ [1 - ri2k{eRk)] + ri2i(eL) + r23i{fL) - ri2k{eRk) = -0.040 ± 0.015, (37) 

where g^^ = -0.0397 d= 0.0003 and gf^ = -0.5064 ± 0.0001. The values quoted for the g- 
parameters are taken from the 2008 Particle Data Group, who point out that the CHARM II 
results [18] dominate the average values. 
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a) b) 



Figure 3: a) SUSY R-breaking process contributing to ql and gR. b) SUSY process contributing 
to gR. Here gauge invariance requires j = 1,3. 



Including Xi2k{^Rk), Eq. (28), in the 2a joint bounds, we find the corresponding upper bound 
to be: 



lA 



121 



100 GeV^ 



EL 



+ 



|A 



2311 



100 GeV\ 



m 



< 0.130. 



(38) 



To put Eq. (38), a bound on the sum of squares of couplings divided by scaled masses, in the 
context of other bounds, we can use the bound from [17], updated to 2008 data [6]: 



121 



100 GeV 



X A 



m 



231 



100 GeV^ 



< 3.0 X 10" 



(39) 



The bound in Eq. (39) combines the experimental bound on the decay rate for r — > eee with 
its representation in the "double coupling dominance convention" for R-parity violating trilinear 
couplings [17]. The representation of the decay involves the sum of squares of five coupling products, 
and the convention, in this case, serves to place the weakest bound on each product by assuming 
all the others are effectively zero. This example, though not in line with our restriction to flavor- 
conserving processes, allows us to discuss the implications for sfermion masses that follow from ^ 
bounds. Because three unknown masses appear in Eqs. (38) and (39), what one can say about 
the implications of the bounds for the A parameters is limited, even if the individual ^ couplings 
entering the two equations are the same. Concisely put, one can say that whenever the sneutrino 
mass satifies 

n 1 .'^n 

(40) 



V2 X 3.0 X 10-4 

the bound of Eq. (38) is more restrictive than that of Eq. (39), which becomes irrelevant. If, 
instead, Eq. (40) is not satisfied, the above bounds have to be considered together, because the 
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hyperbola described by Eq. (39) will cut through the elliptical region defined by Eq. (38), and part 
of the region allowed by Eq. (38) will be prohibited by Eq. (39). 

Unless we invoke some theoretical prejudice about the relative mass scales, we cannot conclude 
more than that. Only if one of the inequalities includes a lower bound, does the combination of 
bounds lead to a general condition on the masses. We will see an illustration of this situation below, 
when considering the combined bounds on v^,^ and V^^e scattering at la. 

3.2.2 Ve + e^ Ve + e 

Turning to the implications of data on the scattering processes + e ^ + ^ and Ve + e ^ 
z7e + e [19, 20, 21], we must consider both high energy data, Ey » m-e, and low energy data, 
Ey rUe. General, model independent analyses of bounds on non-standard interactions from these 
and related neutrino and electron data have recently been carried out for both non-universal and 
flavor-changing new physics interactions [22, 23, 24]. We focus here on the bounds on Iji trilinear 
coupling parameters provided by flavor diagonal elastic V(,e accelerator data at tens of MeV [20] 
and elastic Vet reactor data at several MeV [19]. 

The LSND Collaboration provides a measurement of the total cross section for clastic scattering 
of the electron neutrinos off electrons. Assuming that the final state neutrinos are also electron-type, 
we can use their reported value and the general expression for left handed neutrinos scattering of 
unpolarized electrons to set a limit. The general expression for a,(y — A) ^[ql^Y — A) ^ guiV +Ay\ 
four-fermion interaction differential cross section reads 



da _ 2Gpme 



dT Tx 
and for the total cross section 



(41) 



2meEyGp / 2 , 1 2 1 "^-e 
= \9l+ o9r- 7:^9LgR 



= (10.1 ± 1.5) X 10~'^^EUGcY) cm^, (42) 

where E^ is the neutrino energy in the rest frame of the target electron, and T is the kinetic energy 
of the recoil electron^. The i'ge cross sections follow by interchanging and gR in Eqs. (41) 
and (42). When E^ S> mg, as in the case of the LSND experiment, me/Ey is ignorable, and the 
expression for the cross section simplifies to the familiar high energy form. Including the ^ trilinear 
parameters in the expressions for the gz and gR coupling coefficients, we find 

9T = Q+^w^) [l-rm(eRfc)] (43) 

9r'' = [1 - rukieRk)] + ruiifj-L) + risiifi), (44) 



^In view of the low precision of the experimental uncertainties, we have not included radiative corrections. For a 
study discussing radiative corrections and future possibilities for precision measurements, see [25]. 
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where we have considered a SUSY process hke the one depicted in Figure 3b, in which 1/2 = has 
to be replaced by ui = Ug, X2ji Mji, and j = 2,3, while the Ai2fe-dependence is given by the 
correction to Gp, Eq. (20). In our study of the bound on ri2k{eRk) that follows from the precision 
measurement of muon decay and the renormalized expression for the muon decay formula, we found 
the bounds of Eqs. (25) and (28) respectively at the la and 2a C.L. The corresponding values of 
ri2k{^Rk) are so small that it can be dropped from further discussion. The coupling coefficient 
then has its SM value, and gn is modified from the SM value by the terms that depend on A121 
and A131. Referring to Eq. (42) and (44), we find the bound on the region of trilinear couplings we 
are after: 




121 



100 GeV\^ , 100 GeV\^ 



+ |Ai3i| < 0.66, (45) 



at 2a. 

Before discussing the tie-in of Eq. (45) with other limits, we look next at the independent limits 
set by the results for I7g + e ^ I7g + e from reactor data. In this case, the electron mass-dependent 
terms are important and must be kept. The cross section expression in Eq. (42) is modified by 
interchange of gL and gu for application to the VgC case. 

3.2.3 i7g -f e — > z7e -|- e 

The highest statistics experiment Ve + e ^ Vg + e is still that of Reines, Gurr and Sobel [19]. 
The results are presented as dimensionless factors times the SM charged current, V — A expression 
ay -A, for the cross section for each of two kinetic energy bins. The cross section for a given recoil 
electron energy range is evidently the result of folding the differential cross section with respect to 
electron energy with the Ve flux [26], integrating over neutrino energies and then integrating over 
the recoil kinetic energy. The experimental cross sections reported are cJexp = (0.87 ± 0.25)av-A 
(1.5 MeV < T < 3.0 MeV) and (Jexp = (1.70±0.44)(7y-A (3.0 MeV < T < 4.5 MeV). The quantity 
thus calculated is a function of the ^ parameters, which enter through the coupling g'^^, Eq. (44). 

The theoretical expressions for electron-neutrino and antineutrino scattering involve the same ^ 
couplings. Thus, in the spirit of our multi-parameter, multi-experiment approach, we can combine 
data from LSND [20] and Irvine [19] results for both AT bins in a way similar to what we did for 
Rrij, and Rr in Section 3.1. The resulting constraint at the 2a level is: 



,100 GeV^^ ' .o^r._..x2 




121 



m 



VL 



Y ,100GeV\^ 



At first glance it is surprising that the Ue data, with larger uncertainty, produces tighter constraints 
than the fg data. The source of the added resolving power is the gLQR term in the cross section 
expressions, which plays a significant role in the low energy analysis and increases the sensitivity 
to the variation with respect to the ^ parameters. 
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Though the bound Eq. (46) is consistent with zero at 2(7, at the la level, given current values 
for ql and gn, it is not. This in itself is not of special significance, but it affords the opportunity 
to illustrate added implications when "new ^ -physics" is needed to fill a gap between SM and 
experiment. The joint bound from LSND and Irvine at la yields: 



at la {1.65a), and similarly for Ai3i(fL). 

Now we have the interesting situation that, taking the bound Eq. (47) at face value, we can ask 
under what conditions are they consistent with stringent bounds on related parameters, coming from 
more recent data. Examples we have already studied are the bound |Ai2i|(100 GeV)/me^ < 0.037, 
Eq. (28), and the bound |Ai3i|(100 GeV)/me^ < 0.071, Eq. (31). We assumed here a variant 
of the SCD assumption, corresponding to the largest possible range for individual parameters, 
|Ai2jt| |Ai2i| and |Ai3fc| \^i3i\- 

Prom Eq. (47), scaling by mg^, we can now map out regions of (m^^,mf^) space where the 
constraints from Eqs. (47), (28) and (31) are all satisfied, ignoring for present purposes the mixing 
of la and 2a constraints. The result of this analysis is displayed in Fig. 4. To satisfy the lower 
bound shown in Eq. (47), we see that at least one of the masses must be less than the mass mg^^. 
The smaller the mass becomes, the larger the other must be to satisfy the inequalities. This mass 
information can only be obtained if the strict SCD approach is relaxed, as we have done here. 

Alternatively, taking the bounds one at a time, Eq. (48), we find m^^ < 0.185 (0.34) mg^ or 
ruf^ < 0.355 (0.65) mg^ at la (1.65a-). 

The preceding discussion, summarized in Fig. 4, is offered to illustrate the added power that 
multi-parameter analysis provides to probe ^ parameters. Experiments delivering data with high 
statistics at energies of a MeV or so to study T/gC scattering would sharpen the picture, clarifying 
the possible role of ^ SUSY in this sector of neutrino physics. Here we are considering only low 
energy processes, where the four-fermion efi^ective interactions apply, but at high energies the non- 
local effects of the exchanged particle must be included, directly probing the sfermion masses. This 
possibility is afforded by e+e^ — * vV^ results from LEP [28] and, in the future, possibly 100 GeV 
range i/^e v^e and f^e z/g/i scattering experiments such as those proposed by NuSOnG [29]. 

This concludes our exploration of the multi-parameter effects in purely leptonic processes. Next 

we consider some important constraints from semi-leptonic physics. 

*This mass pattern, if it wore borne out by future experimental constraints, would contradict models that predict 
mg^ < mgj^ , as is the case in minimal supergravity, for example [27] . 




(47) 



Projected onto each parameter one finds. 




< 0.32 




(48) 
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Figure 4: The darkened region shows the values of mg^/m^^ and mg^/mf^ allowed by simultaneous 
application of bounds shown in Eqs. (47), (28) and (31). 



4 Semi-leptonic case 

When R-parity violating interactions are taken into account, charge current and neutral current in- 
teraction generally involve more than one coupling at a time, and in some cases these couplings can 
be large and cancel each other. The lesson we take from the leptonic case is that such degeneracies 
can be removed by considering a subset of experiments characterized by the same R-parity cou- 
plings. Then one bounds the couplings by considering the experimental uncertainties on this subset 
altogether. This is even more evident when we analyze processes that involve the semi-leptonic 
couplings X-j^ of Eq. (11). Contrary to the leptonic and hadronic cases, the couplings X-j^ are not 
required by gauge invariance to have any symmetry in their indices. As a consequence, the number 
of effective couplings entering the Lagrangian is much greater than those appearing in Eq. (13), as 
we have mentioned in Section 2. There are thus more processes that must be used simultaneously 
to bound the couplings. What this also means is that, due to the amazing overall accuracy of the 
SM predictions and the great number of tests, there are many more ways to cut down the allowed 
regions of parameter space. As we will see in the following standard examples, when the availability 
of experiments from which we can draw bounds on a particular coupling increases, the bound on 
the coupling tends to approach the one obtained under the SCD. 
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4.1 Universality in pion and tau decay 

In the cases of semi-leptonic couplings, we can obtain behavior similar in nature to the one depicted 
in Fig. 2a. The ratio: 

^ r(r- ^ TT-Vr) ^ p5M lKd + r-3ifc(4fc)P 

would give in the SCD the 2a bounds lAg^^l < 0.092 {dRk) and \\2ik\ < 0.032 (djik). Here, again, 
the uncertainty on the r lifetime is comparable in magnitude to the one on the branching fraction 
to pions, and has to be taken into account. As in the leptonic case, the simultaneous presence of 
both couplings introduces a two-fold degeneracy. Such degeneracy can be removed by considering 
the ratio [3] 

r(7r- ^ e-pe) j^sM /i , 2 



{ 1 + [r'likidm) - r'2ikidRk)] } • (50) 



r(7r ^ II Vu) 

For the purpose of illustrating our multidimensional approach, it is convenient in this case to follow 
the restriction mentioned in [3], so we use Eq. (50) to effectively place two alternative 2a bounds^: 
either 



(51) 



or 



As we have done in Sec. 3.1, we can combine Eqs. (49) and (50) to place 2a bounds on the 
region spanned by Ag^^j. and Ag^^. The result is shown in Fig. 5a. The allowed region, rescaled to 
the masses of the exchanged squarks, can be enclosed in a box of size {|A3]^j.|(di?jt), lAg^^Kdijfe)} < 
{0.098,0.045}. The resulting 2a bound on X'^^^. reads: 

exactly equal to the one obtained by SCD. 

4.2 Unitcirity of the CKM matrix and forward-backward asymmetry 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix elements are experimentally determined by com- 
paring the rates of decays that involve quarks in the initial state to the rate of muon decay. In 
general, nuclear beta decay is used to determine the value of \Vud\^ while the rates for s ulvi 
and h ^ ului m. K and charmless B decay arc used to determine |V^s| and \Vub\- The R-breaking 
processes involved in these decays are shown in Figure 6. 



^Notc here that Eq. (51) impUos a sum over couphngs and exchanged particles. Ncutrinoless double-beta decay 
places a strong independent bound on A'm, which applies to cLr- and WL-exchange equally [30]. 
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Figure 5: a) 2a bound region on A3^^(dij,fc), y^2iki'^Rk) from Rt-,^ and Rj^ combined data, Eqs. (49- 
50). The allowed region can be enclosed in a box of size {| Ag^^l (dj^fe), \^'2ik\i'^Rk)} ^ {0.098,0.045}. 
b) 2(7 bound region on X'^^j^{dRk), cos{A6^^), \'i2ki^Rk) from R^^^ CKM unitarity (Eq. (55), A12A: = 
0) and FB asymmetry combined data, c) 2a bound region on A'^^^^ {dRk), Ai2fc {eRk) and X'^2k i^Rk)- 
\\'i2k\ is bounded by the FB asymmetry while |Ai2fc| by fi decay in the MS scheme. cos(A^^'') = 
— 1. The allowed region can be enclosed in a box of size {\X[^j^\{dRk),\'^i2k\{eRk),\''^'i2k\(^Rk)} — 
{0.047,0.042,0.036}. 




Figure 6: SUSY processes involved in d — > ueue (j = 1), s — > ueu^ {j = 2), and b ueve (j = 3). 

The unitarity constraint can be imposed on the CKM matrix elements, together with the 
effective Lagrangian of Eq. (11) and a similar one, constructed from Eq. (10), involving the 
product of different couplings. One gets [5]: 
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i=l 



1 + n2fc(eRfc 



+ 



+E 



^llfc^l2fc 



+ 



X'* X' 



(54) 



which becomes at leading order in R-parity breaking, 



i=l 



1 - 2ri2fc(eRfc) + 2r'nfc(4fc) + 2 



|a;,j|a;,Jcos(a^ - 

4V2GFm| 



IK. 



ni.ilA'mlcos(AO 



4^/2GFm2- 

dRk 



\Vub\ 



(55) 



where cos(A0^*) = cos{9us + ^i2fe — ^iife) and cos(A0^^) = cos{6ub + Oi3k — ^uk) are the relative 
phases between the CKM matrix elements and the complex R-parity violating couplings. Using 
Eq. (55) we can place bounds on the A' couplings involved by separation between the right- and 
left-hand side. One can substitute the most recent experimental determination of the central values 
of the CKM matrix element on the right, and use the errors on the unitarity bound on the left at 
the desired level of precision. 

In the literature, Eq. (55) is treated in the SCD, with the additional constraint that R-parity 
couplings and CKM matrix elements are treated as real. With these assumptions we find that the 
most recent data [6] result in the following bounds at 2a: |A'j^^^| < 0.027 (c/ijfe) and |Ai2fe| < 0.028 

In Eqs. (54) and (55) the notation of Eq. (16), to express the sum of moduli squared, has been 
used, together with the correction to Gp from the muon lifetime, Eq. (20). As can be seen, the 
full dependence on the CKM and R-parity violating phases is also indicated. We have adopted the 
Wolfenstein parametrization [31] to express the CKM matrix elements. In this parametrization Vus 
is real, while Vub is not. Nonetheless, measurements of the absolute values of the CKM elements give 
|Kb| ~ 0.004, approximately two orders of magnitude smaller than \Vud\ ~ 0.974 and jV^^I ~ 0.226. 
Thus, the behavior of Eq. (55) is almost independent of A'j^g^., as \Vub\ can be neglected. Taking 
into account the fact that Vus is real and |V^6| is tiny, and neglecting for the moment the SUSY 
correction to Gp, Eq. (55) implies 2a bounds on a three dimensional parameter space spanned by 
l-^iifel' I'^'i2fcl cos(^i2fc — 9iik)- -^iifc) can be bounded by vr-decay, Eq. (50). A'j^gfc can be bounded 
by the forward-backward (FB) asymmetry in fermion pair production reactions e~e'^ — > //, which 
we treat in detail in the next subsection. The 2cr-bound region is shown in Fig. 5b. Note that, 
contrary to the other cases in the paper, here the index k has to be common to the three axes in 
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the picture. One can see that, in spite of the fact that the phases are allowed to take on any values, 
the A' parameters are allowed a slightly larger region when cos(A^^'') = —1. We come back to this 
point at the end of this section. 

4.2.1 Forward backward asymmetry 

The forward-backward asymmetry in fermion pair production has been studied at PEP, PETRA, 
TRISTAN, LEP, and SLC. In order to bound \'i2ki'^Rk) we need charm production, e""e"^ cc. The 
SUSY diagram that contributes to this process is depicted in Figure 7a, with u\ ^ U2, X'uf. ^'uk- 
We assume that the right-handed down squark mass is far enough above the Z-pole that we can 
retain our effective Lagrangians, Eqs. (11) and (13), and use the data in [6], dominated by Z-pole 
measurement s. 

The SM expression for the charm FB asymmetry reads [32]: 



3^0 



where 



+l\r\' {{dial? + Wr? + {Slfn? + {fR9lf\ , (57) 
= -QMr){glgl+9y%-9lg%-g%gl) 

+l\rf WlqI? + {fnQ'kf - {slfn? - {fR9lf\ , (58) 

where gL^R are the usual chiral couplings, Qc = 2/3 is the charge of the charm and 

4V2GfM| 



s- M^ + iMzTz Ve^ 

parametrizes the j — Z interference. The R-parity contribution is obtained by the substitution: 



glgl-^glgl-"-'^^. (60) 

So the correction to the SM reads at lowest order, 



ApB = Ap^ 

where 



2Ao Ai 



(61) 



F{r) = Qc^ir) - Vfglgl (62) 

and r has to be calculated at the Z-pole. By using the standard SU{2) x U{1) expressions for 
gL and gR, and adopting the MS scheme value of sin^ 6w for definiteness, one gets the values: 
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gl = -0.2688, gj^ = 0.2312, gl = 0.3459, g% = -0.1541. We obtain the bound at 2a: 



100 GeV I ■ 



As mentioned above, Fig. 5b sliows tliat allowing for tlie A'^^^. and A'^2fc couplings to have op- 
posite complex phases (cosA^?^* = — 1) slightly extends the allowed regions of parameter spaee 
with respect to the SCD. Furthermore, such an extension becomes significant when we also in- 
troduce the leptonic coupling A12A:, bounded by the experimental limits on the muon liftime in 



the MS scheme, Eq. (28). The 2cj-allowed region in A'^^^, Ai2fc and A'^2a- obtained by simoul- 
taneous combination of the data from CKM unitarity, FB asymmetry in charm production and 



muon decay in the MS renormalization scheme is shown in Fig. 5c. It is enclosed in a box of size 
{\Kik\idRk)AM2k\iem),\K2k\idm)} < {0.047,0.042,0.036}, thus allowing roughly factor of two 
extensions of the parameters with respect to the SCD bounds. The combined analysis furnishes a 
new 2(7 bound on A'^^^^^: 

This striking situation, where three parameters are all allowed to be non-zero and larger than their 
SCD values, is obscured when only one parameter at a time is considered, i.e. SCD is assumed 
uniformly. In principle the second row of the CKM matrix could be used in a similar fashion to 
bound |A^ij.| and |A'22fc|: 

^ \V,,f = 1 - 2rMeRk) + 2 ( ^ ^^^^^[^^ ) + 2^r^2.fe) |K.| , (65) 

i=l \ k dftk / * 

where cos(A^^'^) ~ cos(021A; — ^22fc) in the Wolfenstein parametrization. |A2;^;j| (djik) can be bounded 
by pion decay, Eq. (50). The dependence on Ai2fc (eRk) comes from the bounds on universality of 



the Fermi constant in muon decay. We use, again, the MS bound at 2a, Eq. (28). The weakest 
bound consistent with both these constraints is obtained when cos {A9f) = -1 and reads, at 2a, 

A caveat is necessary at this point, in the sense that Eq. (65) is derived for processes involving 
the production of charmed particles in deep inelastic z^^^-nucleon scattering, with the assumption 
of lepton flavor conservation. This is the standard textbook process used for the determination of 
the CKM couplings \Vcd\ and |V^s| [33]. Such a choice is reflected by the i = 2 index of the X'-j^ 
couplings entering Eq. (65). The use of recent PDG2008 data for the uncertainty affecting the 



^As in Footnote 5, A'm is tightly bounded by neutrinoless double-beta decay [30]. 
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unitarity constraint and for the central values of the CKM matrix elements is not fully consistent 
with this idealized picture. The most recent and precise values given in [6] are obtained through 
a weighted average of different processes, some of which involve external particles of the first or 
third lepton generation. It is clear that the robustness of the bound given in Eq. (66) depends 
strongly on the amount and nature of the weighting involved. Because such detailed knowledge 
and extensive analysis in this regard goes beyond the purposes of this paper, we limit ourselves to 
presenting the bound above, recommending caution in its interpretation. As we will see in Section 
4.4, D° decay alone places bounds on the same (sum of) couplings. We consider those bounds more 
robust. 

Finally, Eq. (55) has the nice feature that it involves the phases of the R-breaking couplings. 
In general such phases are associated with CP violating effects. So we can envisage a strategy that 
would combine additional experiments in the CP violating sector with those that can place bounds 
on the moduli of Iji -couplings like the two above, so that a more thorough restriction of parameter 
space takes place. However we did not find in the literature [5] , nor were we able to create a specific 
example that would help us bound the phases of the couplings involved in this case, namely the 
product A']*^;i.A']^2fc (dRk), in terms of CP violating processes. Some asymmetries in fermion pair 
production at leptonic colliders (Z+Z" fjfj') on and above the Z-pole [34] can be expressed in 
term of non trivial combinations of R-breaking phases like ^{^'ijk^'ij'kKjJ^tjJ'yi'^'^ijK'iJ'l'^^ ^^^'^ 
obvious summation over dummy indices. A detailed and comprehensive study of such processes 
would probably shed light on the phenomenological constraints on CP violating phases. Nonethe- 
less, due to the great number of couplings involved, such a study would have to take into account 
a large number of interactions, many of which cannot be treated as "low energy" processes. This 
clearly exceeds the purposes of this paper, requiring an extensive, separate investigation. 

As anticipated above, and shown in Fig 5b, we have tried to constrain the phase difference 
cos(^iifc — 9i2k) by using Eq. (55), where all the absolute values are bounded by some other 
experiments. We have also tried to constrain the phase cos(^?2ifc ^ ^22fe) with the second row, Eq. 
(65). We found no handle to constrain parameters, as any possible values of the phases are allowed 
by CKM unitarity 

4.3 Atomic parity violation 

We can follow the same technique, and use the bounds on X'^^j^ obtained by CKM unitarity and the 
bounds on \i2k obtained in the MS renormalization scheme to place bounds on X'^j^ from atomic 
parity violation (APV). In the SM the Z-exchange between the electrons and atomic nuclei leads to 
parity violating transitions between particular atomic levels. This has been observed for example 
in the 6S 7S transitions of 55^ Cs [35, 36]. The SM contributions are encapsulated in the weak 
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charge which is defined as [6] 

Ql^' = -2 [{A + Z) Cf *^(n) + {2A - Z) Cf ''(d)] , (67) 

where Z is the atomic number, A the atomic mass, and the coefficients Ci{i) are given at tree level 
by 

C!^\n) = -]^ + ^xw and Cf^{d) = ]^-'^-xw (68) 

The corresponding experimental quantities can be expressed in terms of the SM contributions and 
the ]}l processes depicted in Figure 7 [3]: 

Ci{u) = Cf^\u)[l-n2k{eRk)]-r[^k{dRk) (69) 
Ci{d) = Cf''\d)[l-n2k{eRk)]+r[j^{uLj): (70) 

where we have assumed the R-parity correction to the Fermi constant, Eq. (20). The most recent 
determination of the difference 6Qw = Q'w' ~Qw^ cesium can be found in [6] and its expression 
in terms of ^-couplings reads: 

5Qw = -Qw'ri2k{eRk) + 376/11^(4,,) - 422r[^^{uLj). (71) 




Figure 7: SUSY processes involved in Atomic Parity Violation. 



Again, we can first determine the 2a bounds on the semi-leptonic couplings that one can obtain 
by use of the SCD: |A'^^;.| < 0.051 (Jj^fc) and |AyJ < 0.024 {uij). When both semi-leptonic 
couplings are considered, the region of parameter space that is bounded is two dimensional and 
its shape is similar to that of Fig. 2a. The dependence of 5Qw on the leptonic coupling A12A; 
(e_Rfe) due to Gp-correction introduces an additional direction in parameter space, which becomes 
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three dimensional, decay in one of the renormahzation schemes described in Sec. 3 can be 
used to place bounds on A12A;, while pion decay, Eq. (50), can be used to place bounds on A'^^^^. 
The weakest bound is obtained in the MS scheme, Eq. (28). As we have explained extensively 




Figure 8: a) 2a bound region on \'^-^{uLj), ^uki^Rk) and X[^f,{dfik) from APV, /i decay in the 
MS renormahzation scheme and R^^ combined. The allowed region can be enclosed in a box of 
size {\X'iji\{uLj), \Xi2k\{eRk), \X[^k\idRk)} < {0.055,0.043,0.059}. b) 2a bound region on X'^^ki^Rk), 
X'i2k{dRk) and X'^^ki^Rk) from Dq decay, FB asymmetry and Dg decay combined. The allowed region 
can be enclosed in a box of size {\X'22k\{dRk)-, \^'i2k\^dRk)-, 1^32^1 (dm)} < {0.140,0.034,0.359}. 



in Section 4.1, by simultaneously considering these three processes we can delimit a 2a bounded 
region of parameter space, which we present in Fig. 8a. It can be enclosed in a box of size 
{\^'iji\{^Lj)-, |Ai2fc|(e_Rfc), \X'iijJ\{dRk)} < {0.055,0.043,0.059}, thus allowing only marginal extension 
with respect to the SCD for the A'^^^^ and Xi2k parameters, but roughly a factor of two for the A'^^^ 



parameter. The 2a bound on A'^^ we gather from the combined analysis reads^: 



4.4 D decays 



For our last examples, let us now consider D- and Dg-meson decays. We can implement our 
procedure of taking processes that involve one or more of the couplings we have bounded in the 
previous cases, together with others which are at the moment unbounded, and then use the known 
bounds to restrict the boundaries of the allowed multidimensional parameter space to obtain bounds 
on the remaining couplings. Again we use the averages from experimental data as reported in 
PDG2008 and present bounds at the 2a level. 



^See Footnote 5. 
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We consider the following ratios of branching fractions: i?£)+ = B{D'^ K'^ ^i^u^) / B{D'^ — : 
K^e+Ve) and i?£,o = B{D^ R- ii+v^)/B{D^ R-e+Ve)- Their expression in terms of R- 
breaking semileptonic couplings is given by 

2 

tr , I J \ 



Vcs + r22k{dRk] 



ycs + r[2k{dRk) 



(73) 



Rl^ 



1/1.03 is the reduction due to muon phase-space [37, 5]. The bidimcnsional parameter 
space for \/^'i2k\ vs. \y^22k\ presents a degeneracy very similar to the one depicted in Fig. 2a, where 
the maximum values obtained on the axis correspond to simple use of the SCD: \^'22k\ — ^-^^ {dRk)i 
|A'i2fc| < 0.20 {dRk) for D+ decay and |A'22fc| < 0.10 (J/j^), \\\2k\ < 0-21 (J/jfe) for decay 

It is clear that, of the two ratios we have considered, decay places a tighter bound on these 
couplings. On the other hand, we can bound separately ^i2k by means of the FB asymmetry, Eq. 
(61). The bound is tight enough to cut the bidimensional degeneracy almost entirely. We get, as 
a consequence. 



\\'22k\ < 0.090 



m 



(74) 



, 100 GeV / ' 
ten percent stronger than that obtained by SCD. 

Turning to the D~ t' -\-Vji decays for further constraints, we can bound A32J. in the same 
way, starting from the ratio: 

2 

RnXri^) V{D-^T-Vr) 



^cs + r'22kidRk] 



tjSM 



Vcs + r'22kidRk] 



(75) 



where R^^ = 9.76 accounts for the phase-space suppression. One would get, by SCD use, |A22fc| < 
0.27 {dfik) and |A32;j,| < 0.34 {duk). The combined analysis of Dq decay, Eq. (73), and Eqs. 
(75) and (61) yields the 2a region depicted in Fig. 8b, whose margins are given by the box 
{\^22k\idRk)A'^'i2k\idRk),\^S2k\(dRk)} < {0.140,0.034,0.359}, with only slight extension beyond 
the bounds obtained by assuming SCD. This translates to the 2a bound: 



|A32fcl < 0.29 



d-Rk 



100 GeV, 

This bounds is, again, roughly ten percent stronger than that obtained by SCD [38]. 



(76) 



5 Summary and Conclusions 

In this work, we limited our attention to experimental results from a set of standard leptonic and 
semi-leptonic processes and allowed ^ parameters to vary together, constrained by data at the 
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2cr level, to place bounds on their values. We compared the resulting bounds with those obtained 
from the long-standing procedure of allowing only one parameter to be non-zero at a time, which 
has produced a long, useful list of bounds in the literature over the past twenty years or so. Using 
our different approach, we showed that a joint analysis of different experiments involving the same 
subset of couplings can explore regions of parameter space where the bounds are weakened compared 
to the value set by the SCD procedure. More importantly, the 2cr bounds on individual couplings 
obtained by the combined approach are generally different from those obtained by strict SCD. This 
is due to the fact that almost all processes can by expressed in terms of more than one parameter, 
thus introducing correlations between the couplings and degeneracies in the allowed regions of 
parameter space. The combined-experiments approach helps eliminate these degeneracies and at 
the same time maintains the full parameter space structure. These features provide qualitatively 
different information from that available in the literature, whose results are almost exclusively 
limited to isolating parameters and considering them one at a time. New bounds obtained with 
our approach are given in Table 2, where we present a summary of the results described in the 
preceding sections. 

In the case, we found that the requirement that certain trilinear couplings were non-zero, 
combined with simultaneous constraints involving the same couplings but different sfermion masses, 
we could extract hierarchical relationships among these masses. We illustrated this situation in Fig. 
4, where the la allowed area in the space of "mass ratios" is displayed, and in the paragraphs follow- 
ing Eq. (48), where individual la and 1.65(7 mass bounds are shown. To the best of our knowledge, 
this is the first effort, in the context of purely phenomenological bounds on I/i parameters, to find 
constraints among the sfermion masses. 

In conclusion, we can say that, overall, a richer, more complex picture of parameter space 
and, in most cases weaker bounds on ^ parameters result from a multi-parameter, multi-process 
analysis, compared to the analysis of each parameter in isolation. This conclusion is non-trivial 
when, as the case in expressions we consider, parameters enter as sums of squares, suggesting that 
dropping all parameters but one provides the most conservative limit on each. Nonetheless, since 
wc found the allowed ranges of parameters were larger by at most a factor two, we conclude that 
the SCD approach is a reliable order of magnitude estimate of the upper bounds on the individual 
parameters. At the same time, we conclude that fuller analyses, as exemplified here, are needed 
to search for hints that data arc showing R-parity violation in a region of parameter space where 
several parameters arc non-zero. Finally, such analyses are needed to explore mass relations among 
sparticle masses, which requires disentangling couplings and masses by comparisons of theory with 
data. 

Acknowledgements: We thank Danny Marfatia for many helpful conversations on this work, and 
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A (scale) 


Experiment 


Bound(2CT) 


Corr. A 


SCD Bound 






0.037 


none 


NA 


Ai2i(rriAL) 


Ve{Ve)e 


0.36 


Al3l("lfi) 


0.33 


Al2l(»n-ei) 




0.118 


A23l("^fi) 


0.138 


Al3fe(rneflJ 


Rr 


0.071 


Am(meflJ, A23fc(meBj 


0.048 


Ai3i (rrifi ) 




0.36 


Al2l(W',7z,) 


0.33 


A23fe("lefiJ 


Rj- 


0.066 


Al2fc(m£fiJ, Ai3fc(meBj 


0.051 


A23I ("IfL ) 




0.118 


Al2l(^n-ei,) 


0.138 




C K Munitary 


0.039 


Al2fe(me-flJ,A'i2fc(mj^J 


0.027 




Afb{cc) 


0.027 


none 


NA 




Dq decay 


0.090 




0.10 




ij^ /t^ universal. 


0.040 




0.032 




{t^ /t^ universal. 


0.092 


''^21fe('^eflfc) 


0.092 




APV 


0.045 


Al2fc(me^fc),A'ii;t(mj^J 


0.024 




decay 


0.29 


^22ikK-„J'^mK-„J 


0.34 



Table 2: Summary of constraints on A values with their corresponding mass scale in parenthesis. 
The "Experiment" column gives the measured quantities that are the source of the multi- variable 
bound, the "Bound" column. The "Corr. A" column gives the most directly correlated A determin- 
ing the constraint, while the final column, "SCD bound" gives the value of the bound when all the 
relevant A couplings but the one in the first column are set to zero. The note "none" in a column 
means that only one coupling appears in the relevant expression to compare to experiment. The 
note "NA" means that there is no other coupling to set to zero for the case in this row. 
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